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ABSTRACT 



We present a large sample of fully self-consistent hydrodynamical Nbody/Tree-SPH simulations of isolated dwarf spheroidal galaxies 
(dSphs). It has enabled us to identify the key physical parameters and mechanisms at the origin of the observed variety in the Local 
Group dSph properties. The initial total mass (gas + dark matter) of these galaxies is the main driver of their evolution. Star formation 
(SF) occurs in series of short bursts. In massive systems, the very short intervals between the SF peaks mimic a continuous star 
formation rate, while less massive systems exhibit well separated SF bursts, as identified observationally. The delay between the SF 
events is controlled by the gas cooling time dependence on galaxy mass. The observed global scaling relations, luminosity-mass 
and luminosity-metallicity, are reproduced with low scatter We take advantage of the unprecedentedly large sample size and data 
homogeneity of the ESO Large Programme DART, and add to it a few independent studies, to constrain the star formation history 
of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor and Fornax. For the first time, [Mg/Fe] vs [Fe/H] diagrams derived from 
high-resolution spectroscopy of hundreds of individual stars are confronted with model predictions. We find that the diversity in dSph 
properties may well result from intrinsic evolution. We note, however, that the presence of gas in the final state of our simulations, of 
the order of what is observed in dwarf irregulars, calls for removal by external processes. 

Key words, dwarf spheroidal galaxies - star formation - chemical evolution - galactic evolution 



1. Introduction 

Understanding the dominant physical processes at the origin 
of the dynamical and chemical properties of dwarf spheroidal 
galaxies (dSphs) is challenging. The binding energy of the in- 
terstellar medium of these low mass systems, at the faint end of 
the galaxy luminosity function, is weak. Th e injection of energy , 
due to v iolent explosions of supernovae ("Pekel & Silkl T98q: 
Morietal. 1997; Mac Low & Ferrara 1999; Murakami & Babul 



2005 



1999 : iMori et al.l i2 002; H ensler et al.1 120041: iRicotti & Gnedin 



Kawata et al.. ,2006). or the cosmic UV background 



during reionization (lEfstathioul 119921; iBarkana & Loebl 1 19991: 
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iBullock et ani2000l: iMaver et al.ll2006l) may leave dSphs totally 
devoid of gas and consequently quench their star formation. In 
this picture, the majority of the dSphs are fossils of the reion- 
ization epoch and are ch aracterized by an old stellar population 
(IRicotti & Gnedinll2005h . 

However, observations offer evidence for more complex 
star formation hi stories and reveal a clear v ariety of dwarf 
galactic systems jMateol [19981; lDolphinll2002l) . The spread in 
stellar chemical abundances, and in particular the low [a/Fe] 
values compared to Galactic halo stars at equal metallic- 
ity, are hardly compatible with an early termi nation of star 
forma t ion at the epoch of reionization (e.g.. iHarbeck et all 



20011; IShetroneetalJ 119981 1200 iL 120031; iTolstov et al.l 12003 



Geisler et al.ll2005l: iKoch et alJl2008i examples taken in relation 
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to the galaxies studied in this work). Whilst some dSphs are in- 
deed consi stent wi th rather short star formation episodes, such 
as Sextans dLee et a l. 2003) or Sculptor (^abusiaux et al. 200j|), 
others are characterized by much more extended periods, h ke 
Carina dSmecker-Hane et al1ll996 ; iHurlev-Keller et aT ] ll998h or 
Fornax (IColeman & de Joiig|l2008l) . 

Long durations of the star formation were early advo- 
cated by purely chemical evolution models constrained by 
the dw arf metallicity distribu t ions and color-magnitude dia - 
grams (llkuta & Arimotoll2002l: [Lanfranchi & Matteuccill2004 . 
Subsequent simulations of dwarf galaxies introduced the 
role of the dark mat ter coupled to the stellar feedback 
dFerrara & Tolstovl2000h . and later the full dynamical physics of 
the gas and dar k matter, by means of N-Bod yn-SPH treatment. 
Along this line, iMarcohni etall ( l2006l |2008[) concluded that a 
prolonged (compared to instantaneous) star formation requires 
an external cause for gas removal, which cannot be due to galac- 
tic winds. Intermittent episod es of star formation were at the 
focal point of the analysis by IStinson et al.l ( l2007h . They natu- 
rally arose fr om the alt ernation between feedback and cooling of 
the systems. IValcke et al. (2008) confirmed their self -regulated 
form. These authors also found a grad ual shift of the star fo rma- 
tion towards the inner galactic regions. [kawata et al.l (|2006|) had 
looked for evidence of spatial variation as well, in the form of 
metallicity gradients, but had to stop their simulations at redshift 
1 . Likewise, considering a cosmological box a s initial conditions 
instead of individual halos. lRead et alj ( l2006h stopped their sim- 
ulations early, and focused on the smallest and most metal-poor 
dwarf galaxies. 

In all these works, gas remains at the end of the dSph evo- 
lutions. The resolut ion of this problem constitutes a challenge. 
iMaver et al.l ( l2006l) performed simulations of gas-rich dwarf 
galaxy satellites orbiting within a Milky Way-sized halo and 
studied the combined effects of tides and ram pressure. They 
showed that while tidal stirring produces objects whose stellar 
structure and kinematics resemble that of dSph s, ram-pressure 
strippi ng is needed to entirely remove their gas. [Salvador! et al.l 
( |2008|) proposed a semi-analytical treatment in a hierarchical 
galaxy formation framework and achieved the smallest final gas 
fraction. 

Despite real limitations, such as scarce comparisons with ob- 
servations, incomplete time-evolution, or ad hoc parameteriza- 
tions, we are witnessing a rapid convergence toward understand- 
ing the formation and evolution of dSphs. A critical step forward 
must be undertaken with a large set of simulations to be con- 
fronted with an equally broad sample of data. In particular, the 
chemical imprints resulting from different hypotheses have not 
yet been fully capitalized on. 

The VLT/FLAMES instrument, with fiber links to the 
GIRAFFE and UVES spectrographs, has enabled a revolution in 
spectroscopic studies of resolved stellar populations in nearby 
galaxies. It is now possible to measure the abundances of a 
wealth of chemical elements for more than 100 stars at once. Our 
ESO-Large Programme DART (Dwarf Abundances and Radial 
velocity Team) is dedicated to the measure of abundances and 
velocities for several hundred individual stars in a sample of 
three nearby dSph galaxies: Sculptor, Fornax, and Sextans. We 
have used the VLT/FLAMES facility in the low resolution mode 
to obtain Call triplet metallicity estimates, as we ll as accurate 
radial velocities out to the g alaxies' tidal rad ii ( Tolstoy et al 



as Mg, Ca, O, Ti, Na, Eu. dVenn & Hilll2005l:lLetarte et alj2007l 
Hill et al in preparation; Letarte et al. in preparation). 

In the following, we take advantage of the statistically sig- 
nificant DART sample and data homogeneity, and include some 
recent independent studies, to constrain the star formation his- 
tory of five Milky Way dSphs, Sextans, LeoII, Carina, Sculptor, 
and Fornax. For the first time, populated [a/Fe] vs [Fe/H] dia- 
grams can be confronted with model predictions. Our first goal 
is to establish how well one can reproduce the apparent diver- 
sity of dSph star formation histories in a common scheme. We 
choose to model galaxies in isolation, as this is the only way 
to control the effect of all parameters at play, and to understand 
the dominant physical processes. We will try to see if a com- 
plex star formation history may result from intrinsic evolution 
or if external processes are necessary. We have performed an 
unprecedentedly large number of simulations. Not only do they 
account for the gravity of the dark matter and baryons, but they 
also contain a large number of additional physical mechanisms: 
metal-dependent gas cooling above and below lO'^K, star for- 
mation, SNIa and SNII energy feedback and chemical evolution. 
We focus on the luminosity, star formation history and metallic- 
ity properties of dSphs, rather than on their dynamical proper- 
ties, which turn out to be less constraining. 

The paper is organized as follows: The code and the imple- 
mentation of physical processes are described in Section|2] The 
initial conditions are detailed in Section [3] The presentation of 
the results is split in three parts: Section|4]focuses on the global 
evolution of the galaxies and discusses the main driving param- 
eters, while Section |5] is devoted to the scaling relations. The 
detailed analysis of the chemical properties of Sextans, Leo II, 
Carina, Sculptor, and Fornax are treated in SectionjS] Section |7] 
offers a physical interpretation of the results. Section [8] summa- 
rizes our work. 



2. The code 

We have a d apted the code treeAsph originally developed by 
Serna et alj ( I996I) . with further developments presented in 
Ahmi et all ( l2003l). The chemical evolutio n was intr oduced by 
Poirieret al.l ( 12003. Poirier, PhD thesis) and jPoirier e't al. (200 2|). 
For the sake of simplicity, we recall below the main features and 
the general philosophy of the algorithms. 



radial velocities out to trie galaxies tidal rad ii fl lolstov et al 
2004"; 'Battag ha et all 120061: iHelmi et alj 120061: iBattaglia et aT 



.2008a b). Each of the three galaxies has also been observed at 
high resolution for about 80 stars in their central regions, to ob- 
tain detailed abundances for a range of interesting elements such 



2.1. Dynamics 

All gravitati onal forces are compu ted under the tree algorithm 
proposed by iBarnes &Hutl ( 119861) (see also iHernquistI [19871) . 
This technique is based on a hierarchical subdivision of space 
into cubic cells. One approximates the forces due to a cluster 
of particles contained in a cubic cell and acting on a particle / 
by a quadrupolar expansion of the cluster gravitational poten- 
tial. This is done under the condition that the size of the cubic 
cell is small compared to its distance to the particle /. The ra- 
tio between size and distance must be smalle r than a toleranc e 
parameter, 0, fixed to 0.7 in our simulations (lHernquisllll987l) . 
Indeed, under this condition, the internal distribution of the par- 
ticles within cells can be neglected. Consequently, the number of 
operations needed to compute the gravitational forces between 
particles scales as ~ N log A^, instead of ~ A^^ if one were to 
consider each individual pair of particles. 

The hydrodynamics of the gas is followed in the Lagrangian 
Smooth Particle Hydrodynamics (SPH) scheme. Its allows to 
describe an arbitrarily shaped continuous medium with a finite 
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Fig. 1. Normalized cooling function as a function of temperature 
and metallicity. 



nu mber of part i cles ( Lucy" '1977^, 'Gingol d & Monaghanlll977l) . 
see iMonaghanI (Il992h and PricCj (,2005.) for reviews. Each gas 
particle has its mass spatially smeared out by a smoothing ker- 
nel W (here a spline function). 

Unlike the gravitational forces, which are determined from 
the interactions with all other particles in the system, the hy- 
drodynamical forces result from the contributions of a modest 
number of neighbors. The spatial resolution is determined by 
the smoothing length /i, associated with the particle /, computed 
through the requirement that a sphere of radius /i, centered on 
particle / contains 32 neighbors. 

The integration scheme is the symplectic leapfrog used with 
adaptative time-steps. 

2.2. Cooling 

At temperatures lower than lO ^'K, cooling in prim ordial gas is 
dominated by the molecule j/^. lGalli & Pallal ( ll998h have shown 
that the cooling efficiency of H2 is determined by its mass frac- 
tion xh^ ■ Unfortunately, an accurate computation of xhj is dif- 
ficult and requires to take into account the comple x processes 
of the formation and destruction of H2. Following iMaio et al.l 
(I2007h . we instead fixxH2^o 10 ^. Once the gas is enriched with 
metals, these are important to the cool ing properties. W e con- 
sider oxygen, carbon, silicon and iron (iMaio et al.l l2007l). smce 
they are the most-abundant heavy atoms released during stellar 
evolution, particularly by the SNe II and SNe la, that we follow 
in our simulations. We set the density of the free electrons over 
that of hydrogen, (ne/wH) to 10""* cm"-'. 

Above 10^ K, the cooling function is calculated following 
the m etallicity dependent prescription of [Sutherland & Dopital 
(Il993h . The full normalized cooling function is shown in 
Figure [U for a large range of metallicities. 

A more detailed modeling of the gas cooling is not nec- 
essary, as long as simulations are limited in spatial resolution. 
Indeed, the cooling of the gas is directly dependent on its density. 
Therefore, a limited resolution smoothes the density fluctuations 
of the interstellar medium. 



2.3. Chemical evolution and stellar feedback 

The chemical enrichment of the interstellar medium (ISM) de- 
pends on the interplay between different physical processes. It 
requires us to follow the rate at which stars form, the amount of 
newly synthesized chemical elements, the mass and energy re- 
leased during the different stellar phases, and, finally, the mix- 
ing of the metal-enriched stellar outflows with the ISM. The 
computation of this cycle is done by implementing the original 
equations of chemical evolution formalized by iTinsley, (.1980h . 
as closely as possible. 



2.4. Star formation 

We ad opt the now classical recipe of lKatzl (1 19921) and lKatz et al.l 
(Il996l) . A gas particle becomes eligible for star formation if it is 
collapsing (negative velocity divergence) and //) its density is 
higher than a threshold of psfr = 1.67xl0"^^g/ cm^ . However, we 
do not require the dynamical time to be shorter than the sound 
crossing time (Jeans instability). 

The gas particles, which satisfy the above criteria, form stars 
at a rate expressed by: 



dp* c* 



(1) 



which mimics a Schmidt law (ISchmidllll959h . c* is the dimen- 
sionless star formation parameter, fg is taken as the maximum of 
the local cooling time and the free-fall time. 

For a given time interval Af, a gas particle of mass m„ has a 
probability p* to form a stellar particle of mass m*, where p* is 
defined by: 



Too 

P* = — 



1 - exp|-— Af 



(2) 



The new stellar particle is initially assigned the position and ve- 
locity of its gas progenitor Subsequently, gas and stellar veloc- 
ities are modified in order to conserve both energy and momen- 
tum. 

Each stellar particle represent a cluster of stars, sharing the 
same age and metallicity, whose i nitial mass fun ction (IMF) 
<l>(m) is described by a Salpeter law (ISalpeteij|1955h : 



<E)(m) = 
with X = 



x + I 



m , 



(3) 



-1.35, OT/ = 0.05 Mo andm„ = 50 Mq. 



2.5. Ejecta 

We neglect stellar winds, since they contribute little to the evolu- 
tion of the chemical elements that we consider (magnesium and 
iron), and because the injection power t o the ISM is dominated 
by SN explosions dLeiflierer et al.lll992l) . 

The amount of energy, mass and metals ejected by a stellar 
particle during a time interval At is calculated by considering the 
mass of stars exploding between t and t + At. The dependency of 
the stellar lifetimes on meta llicity is taken into account following 
iKodama & Arimotol(ll997l private communication). 

Hence, the feedback energy released by a stellar particle in 
the time interval [f, f + Af] is: 



A^sN = m-i,fi [nn(f) £11 + «ia(f) -Eia] , 



(4) 



where m* o is the initial mass of the stellar particle, and «n(0 
and nia(f) are the corresponding number of supernovae SNe II 
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and SNe la per unit mass during Af . The energy released by both 
SNe II (£11), and SNe la (EiJ, is set to 10^' erg. 

With mii \ = 8 M0 and »zu u - 50 M© being the lowest and 
highest masses of stars exploding as SNe II, and m{t) being the 
mass of stars with lifetime t, we can express nu(r) as: 



nii(f) - I dm. 

J msix[m{t+At).mii]] 



m 



(5) 



To ca lculate «ia(f). we adopt the model of iKobavashi et alj 
(120001) in which the progenitors of SNe la have main-sequence 
masses in the range of Mp i = 3 Mq to Mp ^ - 8 M©, and evolve 
into C+O white dwarfs (WDs). These white dwarfs can form 
two different types of binary systems (here labeled 0, either with 
main sequence stars or with red giants. Hence: 



«la(0 



r ^ ^ dm\ V' bj f — ^ dm, (6) 



with ^d(>n) the distribution mass function of the companion 
stars, mi l = max[OT(f + Af),M,;d,i], mj^ - min[OT(f), M/,d,u]. The 
lifetime of SNe la is ~ 2-20 and 0.5 - 1 .5 G yr, depending 
on the binary system. See IKobavashi et al] ( l2000h for the lowest 
M, and higher mass M,;d,u of the WDs progenitors and for the 
values of bi which weights the probability of having one or the 
other binary system. 

The supernova feedback energy is released in the form of 
thermal energy only at the end of each dynamical time-step. This 
procedure avoids the thermal energy to be dissipated instanta- 
neously by the strong coolin g above 10^ K, and mimics the blast 
wave shocks of supernovae ( Stinson et alj|2006h . 

The ejected gas mass fraction due to SNe la is given by: 



Amia(f) = mwD nia(0- 



(7) 



with mwD = 1 -38 M© being the mass of white dwarf. 

The mass of each chemical element k ejected by a stellar 
particle is: 



(8) 



where: 

Amii^ 

and 



r-m(t) 

ii(0 = 

*Jma.x[m 

f 



Pk,u(^) 0(m) dm (9) 

max[m(r+A/),mn.i] 
•w(t) 

+ Zk \ (1 - (^(m) - Pk,ii{m)) ^{m)dm 

max[m(A+ A/),/)2ii I ] 



Amkx-M = Amia(f) Pk,ia(m). 



(10) 



(L>{m) is the remnant mass fraction, which is the mass fraction 
of a black hole, neutron star or a white dwarf, depending on the 
initial mass m. Values of a»(m) are taken from (IKobavashi et alj 
I2OOOI) . Zk is the original stellar abundance of the element k, 
Pk,ii{ni) and pk,iJjn) are the stellar yields, i.e. the mass fractions 
of newly produced and ejected element k, com ing from SNe II 
and SN e I, respectively. They are taken from iTsuiimoto et alj 
([T995h . 

Since the stellar particles correspond to star clusters,^e use 
the single st ellar population mass-to-light ratios of iMarastonl 
(Il998ll2005h to calculate their luminosities in V-band. The ef- 
fects of metallicity and age are taken into account. 



3. Initial Conditions 

3. 1 . Mass distribution 

We consider dSphs in isolation. Gas and dark matter are initially 
represented by pseudo-isothermal spheres: 

Pc 



p(r) = 



(11) 



where r is the radius, is the scale length of the mass distribu- 
tion, and Pc the central mass density. The models are truncated at 
'"max- The halo and gas mass distributions differ by their central 
density, Pc.haio and p^gas, respectively. 

As the the total mass inside a radius rmax is linearly depen- 
dent on the central density, there is a proportionality relation be- 
tween the fraction of baryonic matter, /b, and the central densi- 
ties: 



M, 



gas 



(12) 



Mgas + Mhalo Pc.gas + Pchalo 

In the following, we will use pctot = Pc.gas + Pchalo- Using /b and 
Pc.tot, one can simply write: 

Pc.gas — /bPc,tot; 

Pchalo - (l-/b) Pctot- (13) 

Similarly to the gas, the dark matter halo evolves under the 
laws of gravity. 

We consider a core in the initial dark model profile. 
Whilst cosmological simulations predict the formation of 
cuspy dark halos (INav arro et al. 1997; Fukushig e & Makinol 



ipy da 

il997i iMoore et al.1 11998; Spring el et al. 2008, and the ref- 
erences therein), our choice is motivated by observational 
evidences found in normal , low brightness and dwarf 



galaxies dBlais-Ouellette et alj 120011: Ide Blok & Bosmal 12002 : 



'Gentile et al.' '2004' 



iSwaters et al.' '2003' 
i2005; de Blok .2005; .de Blok et al.. .2008> 



^200? 



Spekkens et 
Spano et al. ,2008h . 



Measuring the in ner slope of the dSph p rofiles is very challeng- 
ing, nevertheless. iBattaglia et al] ( l2008ah show that the observed 
velocity dispersion profiles of the Sculptor dSph are best fitted 
by a cored dark matter halo. 

3.2. Velocities and temperature 

The initial velocities are obtained by assuming equilibrium, free 
of any rotation. For a spherical distribution, we can assume that 
the velocity dispersion is isotropic. It can be derived from the 
the second mornent o f the Jeans equation dBinnev & Tremain3 
ll987l:lHernquisl|[T993l) . In spherical coordinates, one writes: 



1 r°° 

crHr)^— dr'p(r')d,,<l>(r'). 
p{r) Jr 



(14) 



The halo velocities are randomly generated in order to fit the 
velocity dispersion cr- at any given radius. 

The temperature T of the gas is deduced from the virial equa- 
tion: 

GM{r) 3k^T 



lim 



P'«h 



(15) 



leading, for a pseudo-isothermal sphere, to: 

r = — — Gpctot'-c, (16) 
3 k-B, 

where k^ is the Boltzmann constant, m^ the proton mass and yu 
the mean molecular weight of the gas. 
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3.3. Initial Parameters 

All simulations start with an initial radius of r^ax - 8 kpc, dis- 
tance at which the gas and dark matter densities are about 1 / 1000 
of the central ones. We consider two different core radii, r^, of 
0.5 and 1 kpc. The choice of the central total density pc.tot (dark 
matter + gas) uniquely determines the initial total mass of the 
system, M,, which we vary over a range of 2 x 10^ to 9 x 10^ Mq. 
We investigate the effect of the initial baryonic mass fraction, 
/b, by varying it from 0.1 to 0.2. Indeed, this helps in disentan- 
gling the influence of the total gravitational potential from that 
of the gas mass. The masses of the gas and halo particles remain 
constant at 1 .4 x 10"* Mq and 9.3x10"* Mq, respectively. The cor- 
responding gravitational softening lengths are 0.1 and 0.25 kpc. 
As a consequence, the simulations start with 4000 to 20'000 par- 
ticles. The variation in number of particles is therefore at most 
a factor 5, hence a factor of 1.7 in spatial resolution, which jus- 
tifies the choice of fixed softening lengths. The star formation 
parameter c* is varied from 0.01 to 0.3. The initial mass of the 
stellar particles is 980 Mq, corresponding to about one tenth of 
the initial mass of the gas particles. 

4. Models 

We performed 166 simulations to understand the role of each 
parameter at play, and to identify a series of generic models re- 
producing the observations. The complete list of the simulations 
is given in Tab. lA.ll IA.2I and IA.3I of Appendix [A] The models 
have been run for 14 Gyr. 

4.1. Modes 

As mentioned in Section [3] the initial sphere of DM+gas is in 
equilibrium under adiabatic conditions. At the onset of the sim- 
ulations, the energy loss due to cooling causes the gas to sinks 
in the potential well and contract. The total potential is deepen 
not only due to the central increase in gas density, but also as a 
consequence of the halo adiabatic contraction. Despite the large 
increase in density, the gas temperature is kept nearly constant 
due to the strong hydrogen recombination cooling above lO'^K 
(see Fig. [T]i- Therefore, the gravitational energy recovered from 
the deepening of the potential well is dissipated nearly instan- 
taneously. For densities above psfr, the evolution of the model 
depends on supernova heating, directly linked to the star forma- 
tion rate. 

Besides this general description common to all simulations, 
we identify three different major regimes. We refer to them 
as "full gas consumption", "outflow" and "self-regulation". For 
each of them. Fig. |2] presents the evolution with time of the star 
formation rate (SFR), the central gas density, the central gravi- 
tational potential and the mass of the gas within 3 kpc from the 
galaxy center. 

4.1.1. Full gas consumption 

In cases where c* is low for a given M-,, the energy released 
by the supernova explosions is unable to counterbalance the 
radiative cooling. As a consequence, the gas keeps on sink- 
ing in the galaxy inner regions and reaches very high densities. 
Stars are formed continuously and at high rate. The model #647 
(c* - 0.05 and M; - 6.6 x IO^Mq) in Figure [2] provides a clear 
example of this regime: steep rise in star formation rate and cen- 
tral gas density. The chemical enrichment of the resulting sys- 
tems is rapid, and their metallicities quickly exceed the highest 




12 3 1 

t [Gyr] 



Fig. 2. Time evolution of the star formation rate, SFR, the gas 
central density, p^, the galaxy central gravitational potential, Oc, 
and the gas mass, Mg within 3 kpc. The black line illustrates the 
case of full gas consumption. The red line stands for the case 
of an outflow. Finally, the blue line indicates the evolution of 
a self-regulated model. The dotted black line delimits the star 
formation threshold density. 



ones measured in dSphs. These models were not investigated 
further 

4.1.2. Outflow 

Stars can be formed at slightly lower densities by increasing c* 
at a given initial mass, or by decreasing the initial mass at fixed 
c*. This is sufficient to stop the drastic accumulation of gas at 
the center Nevertheless, the gas density is still high, and star 
formation is very efficient. When SNe explode, a huge amount 
of energy is deposited in the gas, which in turn is expelled from 
the galaxie's central regions. In parallel, the central potential in- 
creases (it is negative), primarily due to the ejection of the the 
gas, but also due to the ensuing DM halo expansion. The final 
consequence is a strong outflow. A large fraction of the total 
gas mass is ejected beyond a radius of 3 kpc, chosen to be large 
enough compared to the stellar extent of the systems. For clarity, 
we illustrate this regime with M#519 in Fig|2] in which the out- 
flow occurs early in the galaxy evolution: after ~ 2 Gyr, there is 
virtually no gas left. 

4.1.3. Self-regulation 

Dwarf galaxies are formed in a regime of self-regulation, 
characterized by successive periods of cooling and feedback. 
Such intermittent star formation episodes occurring sponta- 
neously in hydrodyn ami cal simulations have been mentioned by 
IStinson et all (|2007|) and lValcke et"aLl (|2008|) . 

M#576 in Fig|2] offers the example of an intermediate mass 
self-regulated system (Mi = 4.39 x 10^ Mq). As usual, the 
first contraction of the gas leads to a peak in star formation 



6 



Y. Revaz et al.: The Dynamical and Chemical Evolution of Dwarf Spheroidal Galaxies 



{t = 0.6 Gyr). The gas expelled by the supernova feedback is di- 
luted at densities below psfr, and the star formation stops. As the 
gas particles cool, they become eligible to star formation again, 
forming a new burst. Star formation occurs at high frequency in 
M#576 (periods between 100 and 200 Myr). It produces a flat 
distribution of stellar ages, mimicking a nearly continuous star 
formation rate (see Fig. [3j. We will show later that the corre- 
sponding chemical signatures are also very homogeneous. 

Contrary to what has been observed bv lStinson et al.l (|2007|) . 
the fluctuation of the SF is not strictly periodic. However, we 
confirm the influence of the total mass M\ on the duration of 
the quiescent periods jValcke et alj|2008l) . Lower mass systems 
(Mi ^ 3 X IO^Mq) are generally characterized by star formation 
episodes separated by longer intervals, up to a few Gyrs. These 
systems exhibit inhomogeneous stellar populations. 



4.2. Driving parameters 

The description of the different regimes of star formation his- 
tories already points out the importance of both the initial total 
mass. Mi, and the star formation parameter, c*. Figure[3]presents 
the stellar age histograms for models with /b = 0.15 and rc = 1. 
From bottom to top, increases by a factor 6. From left to right. 
Mi increases by a factor 3. The highest mass systems are charac- 
terized by a strong predominance of the old stellar populations. 
Decreasing the initial total mass extends the period of star for- 
mation, passing progressively from a continuous to a discrete 
distribution of stellar ages. The role of c* appears secondary, 
distributing slightly differently the difi'erent peaks of star forma- 
tion (position and strength). 

Fig. |4] summarizes the 166 simulations in a diagram of 
and Mi, for different core radii rc and baryonic fractions /b. 
Colors code the final galaxy stellar metallicity ([Fe/H]), com- 
puted as the median of the distribution, since it best traces the 
position of the metallicity peaks in the observations. The size 
of the circles is proportional to the final stellar luminosity in 
the V-band, Ly The small black triangles indicate the simula- 
tions that lead to full gas consumption and have been stopped. 
As described earlier, the latters result from a too small c*. It 
can be avoided, for our purpose, by increasing or decreas- 
ing Mi. Self-regulated systems with limited outflow are found 
for smaller Mi. These tendencies do not dependent on /b and r^, 
which can only slightly modify the interval of mass in which 
a particular regime is valid. For a given c*, ([Fe/H]) increases 
with Mi and similarly, for a given Mi, ([Fe/H]) increases with 
c*. At very low mass, however, ([Fe/H]) is only weakly influ- 
enced by c*. On the contrary, the larger the mass, the smaller 
increase is needed to raise ([Fe/H]). 

The left and middle panels of Fig.|5]display the final galaxy 
stellar metallicity and stellar V-luminosity, respectively, as a 
funct i on of Mj. The Lo cal Group dSphs luminosities (Mateo 
II998l:lGrebel et al.ll2003l) and mean metalUcities (DART) are in- 
dicated with horizontal dotted red lines. The most outstanding 
result is that changing Mi by a factor 4 translates to a change 
in Lv by a factor 100. ([Fe/H]) is varied by a factor ~ 3. By 
comparison, the influence of c* on the galaxy properties appears 
small. In any case, increasing will also help increasing both 
Lv and ([Fe/H]). The consequence of varying Mi is not linear At 
low initial mass, a small increase in mass is sufficient to strongly 
increase Ly and ([Fe/H]), while at larger initial mass, the rela- 
tions saturate and a larger step in mass is necessary. The mass- 
luminosity and metallicity-luminosity relations will be discussed 
in the next section. 



The right panel of Fig. |5] presents the relation between the 
galaxy's mean stellar age and Mi. The influence of looks more 
linear than previously on Lv and ([Fe/H]). As a matter of fact, 
we have seen in Fig. [3] that it plays a role in the stellar age dis- 
tribution. At a given initial mass, determines the length of the 
star formation periods as well as the interval between them. 

As a conclusion, the above analyses stress the primordial im- 
pact of the initial total mass of the systems. Moreover, one can 
clearly identify the range of possible Mi leading to the forma- 
tion of the Local Group dSphs as we observed them today. This 
range is narrow, e.g, a factor 2 centered on Mi > 5 x IO'^Mq. 



5. Global relations 

dSph galaxies follow luminosity-mass and luminosity- 
metallicity relations that are considered as cornerstones to 
understanding flieir f ormation and evolution ( Mated 



1998 



2008 



IWilkinson et al."2006'; 'Gilmore eXaL||2007t IStrigari et al.] 
Geha et al. 2008; Kirby et al. 200§7^ 

In the following discussion, we calculate all physical quan- 
tities (luminosities, masses, abundances) within the radius /?l 
defined as the radius containing 90% of a galaxy's total luminos- 
ity. This choice is guided by the wish to reproduce as closely as 
possible the observational conditions under which these quanti- 
ties are measured. The classical dSphs (as opposed to newly dis- 
covered faint ones) surrounding the MiUcy Way have tidal radi i 
in the range ~ 0.5 kpc to 3 kpc (llrwin & Hatzidimitriou|[T995l) . 
Fixing a constant small aperture for all dSphs would underes- 
timate both light and mass of the largest systems. Since dark 
matter does not necessarily follow light, this would also bias 
the results. The observational estimates of the dSph total masses 
are based on stellar velocity dispersions measured at galacto- 
centric radii as large as possible, thereby directly linked to the 
limits of the visible matter. Although the farthest measurements 
do not always reach the galaxies' tidal radii, their location is 
determined by severe drops in stellar density, ensuring that the 
bulk of the galaxies' light is enclosed. Consequently, we com- 
pare our models to the mass es derived at the outermost velocity 
dispersion profile point (e .g, IWaUcer et al.ll2007l: iBattaglia et al.l 
l2008atlKievna et al.ll2004l) . Ursa Minor is flie only exception to 
this rule. Its mass has been derived from its central velocity dis- 
persion (Mateo 1998). 

Figures |6] and |7] display the relations of the galaxies' mass- 
to-light ratios M/Lv and the median of the stellar metallicity 
distributions ([Fe/H]), together with the total luminosity of the 
model galaxies. The observations are represented in red, with 
squares for the Milky Way satellites and crosses for t he oth- 
ers. In general, the values of ( [Fe/H]) are t aken f rom iMateol 
(1 19981) when available or from Icrebel et all (|2003|) otherwise. 
The mean metallicities of Carina, Fornax, Sculptor, and Sextans 
are calculate d from their m etallicity distributions (iHekni et al.l 
2006; BaUaglia et al.ll200"8al i2006).The mean metal licity of Leo 
II is d erived from the metallicity distrib ution of Bosler et al.l 
(|2007|) . The luminosities are taken from iGrebel et al.l (l2003h . 
with the exception of Draco (iMartin et all 2008 ). The masses 
of Carina, Fornax, Draco, Leo I and Leo II are computed by 
IWalker et al.l ( |2007|) inside rmax- The mass of Sextans corre- 
sponds to the upper limit of lKIevna et al.l(l2004l) . while the mass 
of Ursa M inor comes f rom'Mateoj (Il998l) . Sculptor's M/Lv is 
taken from lBattagUa et al . (2008^^ 

Both M/ Lv and ([Fe /H] ) show very clear log-linear relations 
with Lv! 



logio (M/L) = -0.791ogio(Lv/Lo) + 1.85, 



(17) 
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Fig.3. The distributions of stellar ages between and 14Gyr for a sub-sample of simulations with /b - 0.15 and = 1. 



and 

([Fe/H]) = 0.68 log,,, (U/Lo) - 1.8. (18) 

Despite differences in c^,, and fi,, all simulations nicely 
reproduce the observations, with a reasonably small scatter. This 
stresses once more that a dSph galaxy's total initial mass drives 
most of its evolution. 

To allow deeper insight into the building up of the M/ Ly - Ly 
relation. Fig. |8] distinguishes between dark matter (DM), stars 
and gas. Colors encode the three different initial baryonic frac- 
tions that we have considered, /b - 0.2 (yellow), /b = 0.15 
(green) and /b =0.1 (blue). 

Very naturally, stellar mass scales with the luminosity. At 
fixed luminosity, the dispersion in stellar mass is of the order 
of 5 X 1O''M0. In fact, a more appropriate way to look at this 
panel is to consider the dispersion in luminosity at fixed stel- 
lar mass, since the dispersion in luminosity is a direct conse- 
quence of various distributions in stellar ages and metallicities. 
This dispersion, of the order of 10* L©. It increases slightly for 
larger masses, for which star formation can last longer, inducing 
a larger number of possible age/metallicity combinations. 

As already discussed in Fig. |5] whilst Ly spans nearly 3 or- 
ders of magnitude, the mass of dark matter varies little. This vari- 
ation is much less than one order of magnitude inside /?l and is 
mostly due to the dispersion among the models. A common mass 
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Fig. 4. The effect of variation of the model initial parameters. The full sample of simulations is shown. Each circle or triangle stands 
for one model. The circle size is proportional to the final galaxy luminosity (see the top of the Figure), while colors code the final 
stellar metallicity (median of the metallicity distribution function), following the scale given at the right hand side of the diagrams. 
Each panel encompasses simulations of identical baryonic fraction and initial core radius. The black triangles identify cases of full 
gas consumption. 



scale (inside Ri), around 1-5x10^ Mq, for the dSph total masses 
seems also favored by the observations, although an exact value 
for this limit is difficult to ascertain , given the large uncertain- 
ties of the mass estimates in general (lMateoll998l;lGilmore et alj 
l2007tlStrigarretani2008l) . 

Interestingly, one can now witness the effect of varying the 
initial baryonic fraction. Galaxies with Ly < 3 x IO^Lq exhibit 
identical DM halo masses, whatever fi,, while for larger lumi- 
nosities, the models of lowest /b require larger halo masses in 
order to generate a similar quantity of stars. For fi, = 0.20, the 
DM halo mass is constant over the whole luminosity range. This 
demonstrates that while dark matter plays a crucial role in con- 
fining the gas, the amount of the latter is also important for the 
most massive dSphs; it must reach a critical amount to enable 
their formation. 

The third panel of Fig.[8]displays the mass of residual gas af- 
ter 14Gyr. It constitutes a very small fraction of the total mass, 
and is therefore undetectable at the level of the scaling relations. 
Quite interestingly, its amount is of the or der of the HI mass 
observed in dwarf irregular galaxies (dis) jGrebel et alj|2003h 
and is weakly dependent on the total luminosity. Nevertheless, 



there is a non-intuitive tendency for the most luminous galaxies 
(> 5 X 10^ Lq) to exhibit less gas than the rest of the systems on 
the sequence. However, the most massive galaxies succeed in re- 
taining most of their gas despite the supernova explosions. Less 
than 60% of their gas is ejected, while this fraction lies between 
70% and 90% for less massive systems, in agreement with the 
results o f IValcke et alJ (l2008h . However, the star formation effi- 
ciency is also higher in more massive galaxies. As a result, the 
gas consumption counterbalances the presence of the large gas 
reservoir. We will come back to this in Section|2l 

Since the gas mass is very much constant at around 10^ Mq, 
the smallest galaxies have proportionally more gas than the mas- 
sive ones (see bottom panel of Fig. [8]). Below 3 x 10^ L©, galax- 
ies have more gas left than they have formed stars. Their gas to 
stellar mass ratios reach 100 at the faint end of the dSph model 
sequence. Meanwhile, stars dominate over gas by a factor 5 for 
the most luminous galaxies. 

As a conclusion, the lessons to be taken from Fig. [8] are 
twofold: i) all dSphs are clearly dominated by dark matter 
The lower the mass of the system, the lower the final bary- 
onic fraction, ii) all model galaxies present an excess of gas at 
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Fig. 5. The metallicity, ([Fe/H]), luminosity in V-band, L^, and stellar mean age as a function of the galaxy's total initial mass, after 
14Gyr of evolution. The circles and squares correspond to - 0.5 and = 1, respectively. The colors encode following the 
scale at the right-most panel of the figure. The positions of our four targeted dSphs are indicated with red dotted lines. Plain line 
indicate the mean values, while dashed ones indicate the one sigma deviation. 



the end of their evolution, as ha s been found in s i milar studies 
(Marcolini et al. 2006, 2008; St inson et al.1 120071: IValcke et al.l 
I2OO8I) . As demonstrated earlier, dSphs cannot originate from 
smaller amounts of gas (for a given final luminosity, metallic- 
ity and age). It is necessary to initiate the star formation in the 
observed proportions. However, it is not yet clear how much of 
this gas must be kept in the subsequent phases of the galaxy evo- 
lution. It is clear that the excess of gas, observed in models in 
isolation, must in reality have been stripped some time during 
the galaxy evolution. 

The quantity of gas falls in the range of HI mass observed 
in dis, always found further away than dSphs from their par- 
ent galaxies, thereby bringing another piece of evidence for the 
role that interactions might play in its removal. It might also be 
achieved in a hierarchical formation framework, for which the 
gravitational potentials are initially weaker 

As we have just shown, the global scaling relations are repro- 
duced by our model with impressive ease. In turn, this conveys 
the idea that the global scaling relations do not form a very pre- 
cise set of constraints. They do not reflect the diversity in star for- 
mation histories that we have illustrated. In order to understand 
the individual histories of the Local Group dSphs, one definitely 



needs to go one step further and consider their chemical abun- 
dance patterns, as well as the information that color-magnitude 
diagrams provide on the stellar age distributions. 

6. Generic models 

We will now select and discuss a series of generic models re- 
producing the properties of four Milky Way dwarf spheroidals. 
Carina, Fornax, Sculptor and Sextans. The choice of these mod- 
els is based on four observational constraints: 

1) The dSph total luminosity, which scales with the total 
amount of matter involved in the galaxy star formation his- 
tory. 

2) The metallicity distribution which traces the star formation 
efficiency. 

3) The chemical abundance patterns, in particular that of the 
a-elements, that determines the length and efficiency of the 
star formation period together with its homogeneity. We use 
magnesium for the a-elements. [Fe/H] and [Mg/H] are de- 
rived from high resolution spectroscopy in the central re- 
gions of the galaxies (Hill et al in preparation; Letarte et 
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Fig. 7. The median of the galaxies' metallicity distributions, 
([Fe/H]), as a function of their luminosities Ly. Each circle 
stands for a model. Colors encode the three different baryonic 
fractions that we have considered, fi, = 0.2 (yellow), /b = 0.15 
(green) and fi, = 0.1 (blue). The best fit to our data (Eq. [TtT i 
is indicated by the dashed line. The red squares represent the 
Milky Way's brightest dSph satellites, the red crosses represent 
the MSl's ones. 
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Fig. 8. The relation between galaxy mass and V-band luminosity 
for the the dark matter, the stars and the gas. The lower panel 
displays the gas over stars mass ratio. Colors encode the three 
different baryonic fractions that we have considered, fi, = 0.2 
(yellow), /b = 0.15 (green) and/b = 0.1 (blue). 
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Table [T] lists the initial parameters of our generic models. 
Table |2]lists their final properties, namely the stellar luminosity 
in the V-band, Ly, the mass-to-light ratio, M/Ly, the gas mass, 
Mgas, the median metallicity, ([Fe/H]), the fraction of stars with 
[Fe/H] lower than -3 and, finally, the stellar age distribution di- 
vided in 3 bins; younger than 4Gyr, between 4 and 8 Gyr and 
older than 8 Gyr. 

Fig. [m[T2l and[T3]displav the stellar metallicity distribution, 
the [Mg/Fe] vs [Fe/H] diagrams and the stellar age histograms, 
respectively, after 14 Gyr of evolution. Fig. [T4l shows the stellar 
age-metallicity relations. 

We do not try to match exactly all properties of our target 
dSphs. Instead, we select the four most satisfying models in our 
sample of 166. More specifically, we allow a freedom of a fac- 
tor 2 in luminosity and a shift of a few tenths of dex in peak 
[Fe/H]. It is beyond doubt that we could fine-tune Mi, c*, 
and /b to exactly reproduce the observations. However, this is 
beyond our scope, which remains to test the hypothesis of a for- 
mation framework common to all dSphs. The dependence of the 
results on the model parameters is illustrated in Fig. |9]and[T0l 
We selected 6 models at low (Ly s 0.5 x 10^ Mq) and high 
(Lv s 3.7 X 10^ Mq) luminosities. The baryonic fraction is fixed 
each time, and only and Mi vary. One sees that around a given 
fixed core of observed characteristics (e.g, luminosity, metallic- 
ity), we could run the models on a finer grid of parameters in 
order to match the galaxy properties in all details. This means, 
for example, to exactly reflect the stellar age distribution, or the 
spread in abundance ratios. 

6. 1 . Global evolution 

Fig. [T6l and [Tt] show the stellar surface density, the gas surface 
density and the gas temperature between and 14 Gyr for our 
four selected generic models. Not only do these reproduce the 
properties of Carina, Fornax, Sculptor and Sextans individually, 
but they also depict the full spectrum of evolutions seen in our 
models. The size of each panel is 20 x 20kpc. 

All gas particles initially share the same temperatures, corre- 
sponding to the galaxy virial temperature, as given by Eq.[T6] As 
soon as the simulations start, the gas looses energy by radiative 
cooling. Consequently, the gas density increases in the galaxy 
central regions. Red areas in the gas density maps identify den- 
sities larger than psft, i.e., they mark regions where gas particles 
may be eligible for star formation. When this is indeed the case, 
the newly formed stars are traced by their high brightness in the 
stellar density maps. 

The SN explosions induce temperature inhomogeneities. 
Indeed, the heated central gas expands and generates a wave 
propagating outwards. If the SN feedback dominates the cooling, 
the gas is diluted and the red color vanishes from the center of 
the gas density maps: star formation is quenched. Such quiescent 
periods are characterized by a nearly constant and homogeneous 
central temperature ~ 10'* K, the temperature at which the ra- 
diative cooling counterbalances the adiabatic heating. When gas 
has sufficiently cooled, it can condense again, and star formation 
is induced again. These periods of star formation and quiescence 
alternate at low or high frequency, depending on the mass of the 
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Fig. 9. The variation in [Mg/Fe] vs [Fe/H], [Fe/H] distribution 
and stellar age distribution for 3 models of similar luminosities 
(Lv - 0.5 X 10^ Mq) and identical initial baryonic fraction {fi, - 
0.15). The stellar fraction is defined as the number of stars in 
each bin divided by the total number of stars. 
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(Lv = 3.7x10*' Me) 



galaxy. For low-mass systems, the cooling time is of the order 
of several Gyrs, it is much shorter for more massive ones. We 
will now see how each of these cases translate into dSph stellar 
population properties. 
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Fig. 11. The stellar [Fe/H] distributions of the four generic 
models reproducing Carina, Sextans, Sculptor and Fornax. The 
model outputs are shown in black and the observations with red 
dashed lines. The stellar fraction is defined as the number of stars 
in each bin divided by the total number of stars. 
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Fig. 12. The stellar [IVIg/Fe] vs [Fe/H] diagrams for the 
four generic models reproducing Carina, Sextans, Sculptor and 
Fornax. The observations are sketched with black dots. Their 
error bars are shown. Colors encode the fraction of stars as indi- 
cated by the bar at the top. 



6.2. Sextans 

The stellar popul ation of Sexta ns is dominated by stars 
older than lOGyr dLee et al.ll2003h . With its mean metallicity 
<[Fe/H]) = -2.05 and a V-band total luminosity of 0.5 x 10*" Lq, 
it falls exactly on our luminosity-metallicity relation (Fig. |7j, 
and its properties are reproduced by the model #590 (L = 
0.53 x IO^'Lq, ([Fe/H]) = -2.02 ) that experiences an outflow. 
Sextans' evolution is dominated by an early period of star forma- 
tion lasting about 3 Gyr. The gas surface density map of Fig. [T7] 
shows the dense central region at the origin of the star forma- 
tion burst (f = 1.9 Gyr). At that time, a small but bright stellar 
system is already formed. After this period, the gas is expelled 
and diluted. No star can form until the last Gyr, when the gas 
has sufficiently cooled down to fulfill the star formation crite- 
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Fig. 13. The stellar age distributions of the four generic models. 
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divided by the total number of stars. 
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Fig. 14. The age-metallicity relations of the four generic models. 
Each dot corresponds to a stellar particle. 



riod of quiescence between 7 a nd 3 Gyr ( ISmecker-Hane et all 
119961: iHurlev-Keller et al.|[T998h . Our set of simulations reveal 
that episodic bursts of star formation are intrinsic features of the 
self -regulated low initial masses systems (Mi g 3 x 10** M©, see 
Fig. O. Shortly after a star formation episode, supernovae ex- 
plode, gas is heated and expands. Because its density is low, its 
cooling time is of the order of several Gyrs. This series of well 
spaced-out SF peaks translate into dispersion in [Mg/Fe] at fixed 
[Fe/H]. This spread is more accentuated than for Sextans, due 
to the extended star formation and the longer intervals between 
bursts. 

Fig.|2]shows that Carina, with ([Fe/H]) = -1.82 and Ly = 
0.72 X 10'' Mq, falls above the relation [Fe/H]-Lv defined by our 
models. As already mentioned, we give priority to constraints 
provided by the spectroscopic data, and allow some flexibility in 
Lv Under these conditions, model #533 (Mi = 2.7x10** Mq) pro- 
vides a very reasonable fit to Carina's properties. The three major 
SF episodes are reproduced; 45% of the stars have ages between 
9 and 15 Gyr, 29% have 4 to 8 Gyr and 22% are younger than 
3 Gyr. Whilst Carina is somewhat less luminous than Sextans, 
its mean metallicity is higher. This is due to a lower initial mass, 
locking less matter in stars, but a longer period of star formation, 
obtained by a slightly higher star formation parameter, which 
avoids outflows. 



6.4. Sculptor 

The Sculptor dSph has formed stars early, over a period of a few 
Gyrs. No significant intermediate age population has been found, 
hence excluding star formation within the last ~ 5 Gyr. The ev- 
idence for low a-element enhancement agrees with an extended 
star f ormation and self-enrichment over a period o f at least 2 Gyr 



I2003h 



Babusiaux et"aLl 120051: [Shetrone et al] 120031: iTolstov et al.l 



Sculptor's properties are well reproduced by model #630 
(L = 2.86 X 10^ Lo, <[Fe/H]) = -1.82). The outstanding feature 
of Sculptor, as compared to the two previous dSph, is the small 
dispersion in the [Mg/Fe] vs [Fe/H] diagram (Fig. [T2b . This is 
guaranteed by a strong initial star formation, followed by sub- 
sequent episodes of lower, but smoothly decreasing intensities, 
ensuring the chemical homogeneity of the interstellar medium. 



ria. This last episode of star formation is an artifact of the gas 
retained by our model, as discussed in Section |5] It is however 
negligible compared to the bulk of the population and does not 
influence Sextans' properties. 

Fig. [12] shows that the bulk of Sextans model stars are lo- 
cated at [Mg/Fe] = 0.3, with however a noticeable dispersion at 
lower values. The dispersion originates from the uneven intensi- 
ties of the star formation peaks. They create regions with diverse 
levels of chemical enrichment. When the intensity of star for- 
mation rises again after a period of semi-quiescence, the ejecta 
of new SNe 11 are mixed with material of low ff-element abun- 
dance. Refinement of the model would require a larger sample 
of observed stars at high resolution, particularly to estimate the 
statistical significance of the dispersion in [Mg/Fe]. 

6.3. Carina 

Carina seems to have experienced a complex evolution char- 
acterized by episodic bursts of star formation, with at least 
three major episodes at around 15, 7 and 3 Gyr , and a pe- 



6.5. Fornax 

Fornax is obviously the most challenging case of dSphs. It is the 
most luminous of aU (15.5 x 10^ Lq), it is metal-rich (([Fe/H]) = 
-1.07) and has experienced multiple periods of star formation. 
Indeed, in addition to very old stars, Fornax reveals a domi- 
nant intermediate-age popu lation, as well as stars of ~ 3 - 4 Gyr 
dColeman & de Jongll2008h . 

Fornax is nicely reproduced by the self-regulated model 
#575. Stars are formed in a high frequency series of short bursts, 
lasting between 200 and 500 Myr. The amplitude of these bursts 
progressively decreases with time, until ~ 7 Gyr. Thereafter, part 
of the gas that had been expelled by previous SN explosions, has 
sufficiently cooled to restore star formation, at ~ 3 Gyr, similar 
to what is observed. The succession of the short bursts of sim- 
ilar intensities mimics a continuous star formation and leads to 
an efficient and homogeneous chemical enrichment (([Fe/H]) = 
-1.03). 

We reach a luminosity L = 12.6 x 10^ L©, slightly below the 
observation of Fornax. However, just as in the case of Carina, the 
exact value of the V-band luminosity is very sensitive to the ex- 
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Fig. 15. The [Mg/Fe] vs [Fe/H] diagram, the stellar [Fe/H] dis- 
tribution, the age-metallicity relation and the global stellar age 
distribution of the model #587, compared to the observations of 
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1.0 


0.15 
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4.00e-04 
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0.5 


0.15 


0.050 


1.60e-03 


3.64 


Sculptor 


630 


1.0 


0.20 


0.035 


5.00e-04 


4.12 


Fornax 


575 


1.0 


0.15 


0.200 


6.66e-04 


5.48 



Table 1. Initial parameters of the five generic models. 



act amount of stars formed in the last Gyrs and may easily vary. 
Stars form up to the end of the simulation. This is consistent with 
the observational evidence of a small number of lOOMyr stars 
dColeman & de Jong||2008h . 

Figure [14] show the stellar age-metallicity relations of our 
four generic models. In all cases, we are far from a one to 
one correspondence between age and metallicity, although be- 
low [Fe/H]=-2.5 (but more safely below [Fe/H]=-3), stars are 
generally older than 10 Gyr The dispersion in metallicity is in 
the range 0.5 to 1 dex at ages younger than 10 Gyr. 

6.6. Leo II 

While writing up our results, the abunda nce ratios in a sampl e 
of 24 stars in Leo II have been published dShetrone et al.ll2009l) . 
We keep Leo II separated from the other individual case anal- 
yses, since its abundance ratios have been obtained from lower 
resolution spectra and therefore have larger uncertainties. 

Fig. [TS] illustrates how model #587 can reproduce the star 
formation history of Leo II. The observed metallicity distribu - 
tion, derived from ~ 100 stars, is taken fr om lBosler et al.l(l2007l) . 
and [Mg/Fe] from ( IShetrone et alj|2009l) . Leo II has a luminos- 
ity similar to that of Sextans, but unlike Sextans, Leo II sus- 
tained a low but constant star formation ra te, leading to an impor- 
tant intermediate age stellar population (llVIighell & Rich|[l996l: 
iKoch et al.ll2007l) . The reason for a continuous (low) star forma- 
tion in Leo II's model is found in the slightly smaller initial total 
mass for the same star formation parameter as Sextans; it pre- 
vents the gas outflow. Carina and Leo II have about the same 
mean metallicities. They probably have probably experienced a 



very similar evolution, as indeed seen in Table [2] for their mean 
properties. However, the two galaxies indeed differ: i) by the pe- 
riods of quiescence which are longer in Carina, and ii) by the 
extent of the star formation peaks, which are between 2 to 3 
times longer in Leo II, at intermediate ages. Model #587 pro- 
duced a higher final luminosity than observed, due to the pres- 
ence of residual gas and therefore recent star formation, as in 
all the models. It is important to recall that we do not aim at 
reproducing Leo II properties in fine details, but instead check 
whether its main features exist in the series of models we ran. In 
this respect, model #587 demonstrates that a low level of rather 
continuous star formation is indeed achievable. 

The extension of the sample of stars with spectroscopic data 
allowing measurements of abundance ratios would greatly help 
the precise identification of Leo II model. This is particularly 
true at [Fe/H] > -1.6 which would confirm or not the apparent 
plateau in [Mg/Fe]. Even more decisive and more general, one 
needs to observationally confirm or disprove the dispersion in 
[Mg/Fe] at low and moderate metallicities in low mass dSphs, as 
it signals the discrete nature of star formation as advocated by 
contemporary models. 

Table [1] gives the fraction of very-metal poor ([Fe/H]< -3) 
stellar particles in the generic models. The general trend of our 
series of 166 models is that the more extended the star formation 
period, the smaller this fraction. This is clearly illustrated by the 
models of Carina and Fornax, as compared to those of Sextans 
and Fornax. Although the tail of very metal-poor stars is small 
(from~ 1 to 6%), it is still large r than suggested by the observa- 
tions to date jHelmi et al.l2006l) . In order to properly address this 
.particular issue, one needs i) to investigat e th e impact of more 
sophisticated IMFs, such as the one of iKroupa et al. (19^, for 
which the number of low-mass stars is smaller than for a Salpeter 
IMF, and ii) to introduce the peculiar features of star formation 
at zero metallicity. The conditions of transition between popula- 
tion III and population II star formation depends on the galaxy 
halo masses, and so is the galaxy metallicity floor dWise & Abell 
l2008l) . This opens a independent new field of investigation. For 
the time being, we check our model predictions against popula- 
tion II properties. Interestingly, the outflows generated by these 
very first generation of stars offers an alternative solution to get 
rid of the final excess of gas. 



6. 7. Mass profiles and velocity dispersions 

Although our primary goal is not to match the morphology of 
the dSphs, we have considered the shape of the final modeled 
dSph profiles within R^. Fig. [T8]displays the mass density pro- 
files of the stellar and dark components for our four main targets. 
Each model is successfully fitted by a Plummer model, reveal- 
ing a flat inner profile, in agreement with the observations. The 
ratio between the stellar and dark core radii rjt^"/rj!^'° ranges be- 
tween 0.2 and 0.5, indicating a spatial segregation of the star s 
relative to the dark halo, as discussed in dPenarrubia et alJ2008l) . 
Additionally, the line-of-sight stellar velocity dispersion pro- 
files are flat in the enclosed within the core radii, with 6 to 
lOkm/s, corresponding to the least and the most massiv e dSph, 



respec t ively, hence well wi t hin the observed range (Wal ker et al 



2006at iMufloz et al.ll2006t iBattadia et alJl2008at iWalker et al 
2006bh . 



7. Discussion 

We have shown that the total initial mass Mi and are the 
two driving parameters for the evolution of dSphs. Our models 
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26.4 
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0.0 
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83 


0.32 
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0.0 


96.4 
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575 


12.76 


27 


0.12 


-1.03 


1.2 


14.4 


10.6 
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Table 2. The final properties of the five selected generic models. (Ly) : stellar luminosity in V-band, (M/Ly) : mass-to-light ratio, 
(Mgas) : gas mass, (([Fe/H])) : median of the metallicity distribution function, ([Fe/H] < -3.0) : fraction of stars having a metallicity 
lower than -3, (< 4Gyr) : fraction of stars younger than 4Gyr, ([4, 8] Gyr) : fraction of stars having ages between 4 and 8 Gyr, 
(> 8 Gyr) ; fraction of stars older than 8 Gyr. 
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Fig. 18. The stars (blue circles) and dark halo (red triangles) den- 
sity profiles for the four generic models. Each profile is fitted by 
a Plummer model. The positions of the core radii are indicated 
by the dashed lines. Their values are given in the left bottom 
corner of each panel. 



which describe galaxies in isolation, i.e., not yet taking into ac- 
count interactions or accretions that inevitably occur in a ACDM 
Universe, already lead to the observed variety of star formation 
histories. Decreasing the initial total mass of the galaxies sepa- 
rates the periods of star formation. 

The marked differences found between high- and low mass 
systems is a key point in understanding the formation of dSphs. 
In order to illustrate our findings, we now focus on cooling 
time and gas temperature. The top panel of Fig. [19] displays 
the galaxy's virial temperature as a function of its initial mass 
(Eq.[T6). 

The green horizontal band indicates the temperature at which 
the cooling time of the gas is equal to 100 Myr, for two threshold 
densities, psf, and lO^^'psfr that are representative of the range in 
gas density observed in the course of our simulations. The value 
of 100 Myr is chosen to be short compared to the galaxy life- 
times. The lower border of the green band corresponds to psfr 
while the upper ones correspond to lO^^psfr- The intersection 
between the green band and the curves corresponds to a total 
mass ranging between 3 and 4x10** M©, delimiting two regimes 
distinguished by long and short cooling times, as can be mea- 
sured in the bottom panel of Fig. \T9\ The cooling times of psfr 
and lO^-'psfi- as a function of the galaxy initial masses are shown 
in solid and dashed lines, respectively, in the bottom panel of 
the figure. They are calculated by taking the virial temperature 



of the corresponding galaxy mass. The gas contained in galax- 
ies of total masses larger than 4 x 10** M©, with corresponding 
virial temperatures larger than 10"* K is characterized by a short 
cooling time (r <100Myr) at all densities. This short cooling 
time is due to the strong radiative cooling induced by the re- 
combination of hydrogen above 10"* K. A galaxy which is in this 
cooling regime will see its gas loosing a huge amount of energy, 
sinking in its central regions and inducing star formation. This is 
precisely what happens in the model #575 where stars are con- 
tinuously formed, generating a metal-rich Fornax-like system. 
By decreasing the total mass from 4 to 3 x 10*^ M©, the virial 
temperature falls below the peak of hydrogen ionization. As a 
consequence, the cooling function drops and the cooling time is 
increased by nearly three orders of magnitude. Below 3x10** Mq, 
the cooling time is long, since the loss of energy by radiative 
cooling is strongly diminished. In this regime, star formation can 
only occur in episodic bursts, separated by periods correspond- 
ing to the mean cooling time. Carina gives a clear example of 
such a case. For even lower mass systems (Mtot < lO'^M©), the 
cooling time is longer than 10 Gyr, and therefore no stars are 
expected to form. 

From our simulations, a total mass of ~ 4 x 10** M© leads 
to a luminosity of ~ 3 X 10'' L©. Not surprisingly, this lumi- 
nosity corresponds to the critical one found in Section |5] below 
which the galaxies have more gas left than they have formed 
stars. The sharp cutoff of the cooling function below 4x10^ Mq 
explains very nicely why a small decrease in mass induces a 
large drop in luminosity, seen in the middle panel of Fig. [5] A 
direct consequence is the constancy of galaxy masses within /?l 
(~ 1 - 5x 10^ Mq) over the wide range of dSph l uminos i ties, a s 
discussed in Section[5](F ig.[8]i and suggested bv lMated ( Il998l) . 
iGilmore et all (l2007h and lStrigari et alj (l2008h . 

The next step in simulating the formation and evolution of 
dSph galaxies is to study the impact of external processes re- 
sulting from a ACDM complex environment, like tidal mass 
stripping, ram stripping, dark ma tter and gas accretion, a s well 
as reionization and self-shielding. iPefiarrubia et al.l ( |2008l) show 
that tidal mass stripping strongly depends on the spatial dis- 
tribution of stars relative to the dark halo extent. If stars are 
very concentrated in large halos, dSphs will be more resilient 
to tidal disruption. Dwarf spheroidal galaxies may need to loose 
nearly 90% of their mass before their star f ormation and there- 
fore chemical evolution start to be affected. iMaver et al.l (l2006h 
tackle all the above listed questions and show how their mecha- 
nisms can be interlaced. As an example, in their models gravita- 
tional tides help ram pressure stripping by diminishing the over- 
all potential of the dwarf, but tides induce bar formation mak- 
ing subsequent stripping more difficult. Reionization prevent star 
formation, but self-shielding helps cooling. Clearly, understand- 
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Fig. 19. Top: The Virial temperature as a function of the total 
initial mass computed from Eq. [16] with rc = 1 kpc. The lower 
border of the green band gives the Virial temperature of gas with 
density equal to psfi - The upper border indicates the Virial tem- 
perature of gas with density equal to lO^^'pst, . The cooling time 
is of the order of 0. 1 Gyr. Bottom: Cooling time of the gas as a 
function of the total mass, assuming a Virial temperature given 
by Eq.[T6l The case of a gas density equal to psfi is shown by the 
soUd line, while the density of lO^^'psfr is shown by the dashed 
line. 

ing the impact of these physical processes on our specific mod- 
eling requires dedicated simulations. 

8. Summary 

We have investigated the formation and evolution of dSph galax- 
ies, which form a specific class among dwarf systems. Indeed, 
they reach the highest metallicities at fixed luminosity and are 
devoid of gas. 

We performed 166 self-consistent Nbody/Tree-SPH simula- 
tions of systems initially consisting of dark matter and primor- 
dial gas. We have only considered galaxies in isolation. This has 
enabled us to identify the dominant physical ingredients at the 
origin of the observed variety in dSph properties. It has also al- 
lowed us to distinguish which of those are due to the galaxy 
intrinsic evolution and for which interactions are required. The 
diversity of star formation histories of the MiUcy Way dSphs, 
Carina, Leo II, Sextans, Sculptor and Fornax have been success- 
fully reproduced in a single formation scheme. 

• The crucial parameter driving the dSph evolution is the to- 
tal initial mass (gas + dark matter). To a smaller extent the star 
formation parameter, c*, influences the final galaxy properties. 
In particular, it governs the stellar age distribution by modifying 
the time intervals at which star formation occurs. Since there is 
no physical reason for changing from galaxy to galaxy, we un- 
derstand its variation, if needed to reproduce the observations, as 
an indirect evidence for different interaction histories. In a hier- 
archical galaxy formation scheme, our initial masses correspond 
to the halo masses that must be reached along the merger tree 
before ignition of the bulk of the star formation. Besides, the 
chemical abundance ratios constrain the timescale of these halo 
mergers. For example, we have shown that Fornax could not be 
formed without a high initial total mass. This could in principle 
mean an already long accretion history. Nonetheless, its chemi- 



cal homogeneity requests that this should be achieved at a very 
early epoch, before star formation in independent smaller halos 
could have widened its abundance patterns. 

• Star formation occurs in series of short periods (a few hun- 
dreds of Myr long). A high frequency SF mode is analogous 
to a continuous star formation, while a low frequency SF mode 
produces well separated bursts, easily identified observationally 
as in the case of Carina. The period between the star formation 
events is governed by the mass of the systems. It is caused by the 
dependence of the gas cooling time on the galaxy mass; the more 
massive, the shorter it is. This mass dependency is explained 
by the drop of the cooling function below 10"* K. Systems less 
massive than 3 x 10*^ have a virial temperature below 10^ K 
and are characterized by weak cooling. Only episodic periods 
of star formation are expected. On the contrary, massive sys- 
tems (Mtot > 4 X 10'^ Mq) have a virial temperature above lO"* K. 
They are expected to form stars continuously and generate more 
metal-rich dSphs. 

• The dSph scaling relations ([Fe/H], M/L) are reproduced 
and exhibit very low scatter We have shown that there is a con- 
stant final total galaxy mass over the wide range of dSph lumi- 
nosity, which appears as a direct consequence of the cooling time 
-mass relationship. 

• The a-elements were traced by magnesium. The dispersion 
in the [Mg/Fe] vs [Fe/H] diagrams is negligible for continuous 
and/or efficient star formation. It is enhanced when the star for- 
mation occurs in bursts separated by periods of quiescence, or 
when the star formation occurs in peaks of uneven intensities. 
Therefore, it is observed in galaxies as different as Sextans and 
Carina. Our study strengthens the need for large and homoge- 
neous samples of stellar spectra obtained at high resolution, that 
are necessary to constrain and improve the models. They call for 
large observing programs dedicated to the chemical signatures 
at low metallicities in order to establish the level of homogene- 
ity of the interstellar medium in the early phases of the galaxy 
evolution. 

• The fraction of very metal-poor stars in our generic models 
ranges from ~ 1 to 6%. Although this is a rather small amount, 
it is still larger than suggested by the observations. In order to 
properly address this particular issue, one needs to carefully in- 
vestigate the effect of more sophisticated IMFs and the peculiar 
features of star formation at zero metallicity. 

• Gas is found in all our model galaxies after 14 Gyr (~ 
IO'Mq). It clearly points to the need for interactions that would 
strip the gas in the course of the dSph evolution. The possibility 
that gas is gradually expelled as halos merge in a hierarchical 
formation scenario also deserves careful consideration. 
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658 


1.0 


0.20 


0.050 


8.00e-04 


6.59 


1.32 


662 


1.0 


0.20 


0.050 


6.60e-04 


5.48 


1.10 


405 


1.0 


0.20 


0.050 


5.00e-04 


4.12 


0.82 


409 


1.0 


0.20 


0.050 


A r\n r\ a 

4.00e-04 


3.29 


0.66 


563 


1.0 


0.20 


0.050 


3.75e-04 


3.09 


0.62 


555 


1.0 


0.20 


0.050 


3.50e-04 


2.88 


0.58 


454 


1.0 


0.20 


0.050 


3.00e-04 


2.47 


0.49 


624 


1.0 


0.20 


0.065 


5.00e-04 


4.12 


0.82 


620 


1.0 


0.20 


0.075 


5.00e-04 


4.12 


0.82 


667 


1.0 


0.20 


0.100 


1 .OOe-03 


8.24 


1.65 


659 


1.0 


0.20 


0.100 


8.00e-04 


6.59 


1.32 


663 


1.0 


0.20 


0.100 


6.60e-04 


5.48 


1.10 


571 


1.0 


0.20 


0.100 


6.50e-04 


5.35 


1.07 


570 


1.0 


0.20 


0.100 


6.00e-04 


4.94 


0.99 


565 


1.0 


0.20 


0.100 


5.50e-04 


4.53 


0.91 


521 


1.0 


0.20 


0.100 


5.00e-04 


4.12 


0.82 


622 


1.0 


0.20 


0.100 


4.75e-04 


3.91 


0.78 


618 


1.0 


0.20 


0.100 


4.50e-04 


3.71 


0.74 


522 


1.0 


0.20 


0.100 


A r\n n a 

4.00e-04 


3.29 


0.66 


564 


1.0 


0.20 


0.100 


3.75e-04 


3.09 


0.62 


554 


1.0 


0.20 


0.100 


3.50e-04 


2.88 


0.58 


523 


1.0 


0.20 


0.100 


3.00e-04 


2.47 


0.49 


668 


1.0 


0.20 


0.200 


1 .OOe-03 


8.24 


1.65 


660 


1.0 


0.20 


0.200 


8.00e-04 


6.59 


1.32 


664 


1.0 


0.20 


0.200 


6.60e-04 


5.48 


1.10 


562 


1.0 


0.20 


0.200 


5.00e-04 


4.12 


0.82 


559 


1.0 


0.20 


0.200 


A r\n n a 

4.00e-04 


3.29 


0.66 


578 


1.0 


0.20 


0.200 


3.75e-04 


3.09 


0.62 


573 


1.0 


0.20 


0.200 


3.50e-04 


2.88 


0.58 


572 


1.0 


0.20 


0.200 


3.00e-04 


2.47 


0.49 


669 


1.0 


0.20 


0.300 


1 .OOe-03 


8.24 


1.65 


661 


1.0 


0.20 


0.300 


8.00e-04 


6.59 


1.32 


665 


1.0 


0.20 


0.300 


6.60e-04 


5.48 


1.10 


632 


1.0 


0.20 


0.300 


5.00e-04 


4.12 


0.82 


699 


1.0 


0.20 


0.300 


3.00e-04 


2.47 


0.49 


623 


1.0 


0.15 


0.010 


6.60e-04 


5.48 


0.82 


621 


1.0 


0.15 


0.025 


5.30e-04 


4.39 


0.66 


649 


1.0 


0.15 


0.050 


1 .OOe-03 


8.24 


0.82 


647 


1.0 


0.15 


0.050 


8.00e-04 


5.60 


0.99 


585 


1.0 


0.15 


0.050 


6.60e-04 


5.48 


0.82 


586 


1.0 


0.15 


0.050 


5.30e-04 


4.39 


0.66 


641 


1.0 


0.15 


0.050 


4.65e-04 


3.83 


0.57 


587 


1.0 


0.15 


0.050 


A r\r\ f\ A 

4.00e-04 


3.29 


0.49 


648 


1.0 


0.15 


0.050 


2.60e-04 


2.19 


0.33 


650 


1.0 


0.15 


0.100 


1 .OOe-03 


8.24 


0.82 


569 


1.0 


0.15 


0.100 


o r\r\ n a 

8.00e-04 


5.60 


0.99 


566 


1.0 


0.15 


0.100 


6.60e-04 


5.48 


0.82 


567 


1.0 


0.15 


0.100 


5.30e-04 


4.39 


0.66 


568 


1.0 


0.15 


0.100 


A r\n n a 

4.00e-04 


3.29 


0.49 


579 


1.0 


0.15 


0.100 


2.60e-04 


2.19 


0.33 


651 


1.0 


0.15 


0.200 


1 .OOe-03 


8.24 


0.82 


574 


1.0 


0.15 


0.200 


8.00e-04 


5.60 


0.99 


575 


1.0 


0.15 


0.200 


6.60e-04 


5.48 


0.82 


576 


1.0 


0.15 


0.200 


5.30e-04 


4.39 


0.66 


584 


1.0 


0.15 


0.200 


4.50e-04 


3.71 


0.56 


577 


1.0 


0.15 


0.200 


4.00e-04 


3.29 


0.49 


580 


1.0 


0.15 


0.200 


2.60e-04 


2.19 


0.33 


652 


1.0 


0.15 


0.300 


1. OOe-03 


8.24 


1.24 


646 


1.0 


0.15 


0.300 


8.00e-04 


6.59 


0.99 


633 


1.0 


0.15 


0.300 


6.60e-04 


5.48 


0.82 



Table A.l. Model parameters for the complete set of simulations. 
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[KpCJ 








LIU iVlo/KpC J 


[lU IVIqJ 


r 1 A/i 1 
[lu ivioj 


634 


1.0 


0. 


15 


0.300 


5.30e-04 


4.39 


0.66 


644 


1.0 


0. 


15 


0.300 


4.00e-04 


3.29 


0.49 


645 


1.0 


0. 


15 


0.300 


2.60e-04 


2.19 


0.33 


596 


1.0 


0. 


10 


0.050 


l.OOe-03 


8.24 


0.82 


597 


1.0 


0. 


10 


0.050 


8.00e-04 


6.59 


0.66 


598 


1.0 


0. 


10 


0.050 


6.00e-04 


4.94 


0.49 


605 


1.0 


0. 


10 


0.050 


4.00e-04 


3.29 


0.33 


599 


1.0 


0. 


10 


0.100 


l.OOe-03 


8.24 


0.82 


600 


1.0 


0. 


10 


0.100 


8.00e-04 


6.59 


0.66 


619 


1.0 


0. 


10 


0.100 


7.00e-04 


5.76 


0.58 


601 


1.0 


0. 


10 


0.100 


6.00e-04 


4.94 


0.49 


606 


1.0 


0. 


10 


0.100 


4.00e-04 


3.29 


0.33 


602 


1.0 


0. 


10 


0.200 


l.OOe-03 


8.24 


0.82 


603 


1.0 


0. 


10 


0.200 


8.00e-04 


6.59 


0.66 


604 


1.0 


0. 


10 


0.200 


6.00e-04 


4.94 


0.49 


607 


1.0 


0. 


10 


0.200 


4.00e-04 


3.29 


0.33 


670 


1.0 


0. 


10 


0.300 


l.OOe-03 


8.24 


0.82 


635 


1.0 


0. 


10 


0.300 


8.00e-04 


6.59 


0.66 


636 


1.0 


0. 


10 


0.300 


6.00e-04 


4.94 


0.49 


692 


1.0 


0. 


10 


0.300 


4.00e-04 


3.29 


0.33 


528 


0.5 


0.20 


0.010 


2.00e-03 


4.55 


0.91 


529 


0.5 


0.20 


0.010 


1.60e-03 


3.64 


0.73 


530 


0.5 


0.20 


0.010 


1.20e-03 


2.73 


0.55 


676 


0.5 


0.20 


0.050 


4.00e-03 


9.11 


1.82 


673 


0.5 


0.20 


0.050 


3.20e-03 


7.28 


1.46 


672 


0.5 


0.20 


0.050 


2.50e-03 


5.69 


1.14 


525 


0.5 


0.20 


0.050 


2.00e-03 


4.55 


0.91 


626 


0.5 


0.20 


0.050 


1.80e-03 


4.10 


0.82 


526 


0.5 


0.20 


0.050 


1.60e-03 


3.64 


0.73 


527 


0.5 


0.20 


0.050 


1.20e-03 


2.73 


0.55 


677 


0.5 


0.20 


0.100 


4.00e-03 


9.11 


1.82 


671 


0.5 


0.20 


0.100 


2.50e-03 


5.69 


1.14 


531 


0.5 


0.20 


0.100 


2.00e-03 


4.55 


0.91 


532 


0.5 


0.20 


0.100 


1.60e-03 


3.64 


0.73 


533 


0.5 


0.20 


0.100 


1.20e-03 


2.73 


0.55 


627 


0.5 


0.20 


0.200 


4.00e-03 


9.11 


1.82 


628 


0.5 


0.20 


0.200 


3.00e-03 


6.83 


1.37 


629 


0.5 


0.20 


0.200 


2.50e-03 


5.69 


1.14 


608 


0.5 


0.20 


0.200 


2.00e-03 


4.55 


0.91 


609 


0.5 


0.20 


0.200 


1.60e-03 


3.64 


0.73 


610 


0.5 


0.20 


0.200 


1.20e-03 


2.73 


0.55 


678 


0.5 


0.20 


0.300 


4.00e-03 


9.11 


1.82 


675 


0.5 


0.20 


0.300 


3.20e-03 


7.28 


1.46 


640 


0.5 


0.20 


0.300 


2.50e-03 


5.69 


1.14 


637 


0.5 


0.20 


0.300 


2.00e-03 


4.55 


0.91 


691 


0.5 


0.20 


0.300 


1.60e-03 


3.64 


0.73 


690 


0.5 


0.20 


0.300 


1.20e-03 


2.73 


0.55 


625 


0.5 


0. 


15 


0.010 


2.13e-03 


4.85 


0.73 


643 


0.5 


0. 


15 


0.025 


2.13e-03 


4.85 


0.73 


684 


0.5 


0. 


15 


0.050 


4.00e-03 


9.11 


1.37 


680 


0.5 


0. 


15 


0.050 


3.20e-03 


7.28 


1.09 


588 


0.5 


0. 


15 


0.050 


2.66e-03 


6.06 


0.91 


589 


0.5 


0. 


15 


0.050 


2.13e-03 


4.85 


0.73 


590 


0.5 


0. 


15 


0.050 


1.60e-03 


3.64 


0.55 


591 


0.5 


0. 


15 


0.050 


1.20e-03 


2.73 


0.41 


685 


0.5 


0. 


15 


0.100 


4.00e-03 


9.11 


1.37 


681 


0.5 


0. 


15 


0.100 


3.20e-03 


7.28 


1.09 


592 


0.5 


0. 


15 


0.100 


2.66e-03 


6.06 


0.91 


593 


0.5 


0. 


15 


0.100 


2.13e-03 


4.85 


0.73 


594 


0.5 


0. 


15 


0.100 


1.60e-03 


3.64 


0.55 


595 


0.5 


0. 


15 


0.100 


1.20e-03 


2.73 


0.41 


686 


0.5 


0. 


15 


0.200 


4.00e-03 


9.11 


1.37 


682 


0.5 


0. 


15 


0.200 


3.20e-03 


7.28 


1.09 


611 


0.5 


0. 


15 


0.200 


2.66e-03 


6.06 


0.91 


612 


0.5 


0. 


15 


0.200 


2.13e-03 


4.85 


0.73 


613 


0.5 


0. 


15 


0.200 


1.60e-03 


3.64 


0.55 


693 


0.5 


0. 


15 


0.200 


1.20e-03 


2.73 


0.41 


687 


0.5 


0. 


15 


0.300 


4.00e-03 


9.11 


1.37 


683 


0.5 


0. 


15 


0.300 


3.20e-03 


7.28 


1.09 


679 


0.5 


0. 


15 


0.300 


2.66e-03 


6.06 


0.91 
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638 


0.5 


0. 


15 


0.300 


2.13e-03 


4.85 


0.73 


695 


0.5 


0. 


15 


0.300 


1.60e-03 


3.64 


0.55 


694 


0.5 


0. 


15 


0.300 


1.20e-03 


2.73 


0.41 


537 


0.5 


0. 


10 


0.025 


4.00e-03 


9.11 


0.91 


538 


0.5 


0. 


10 


0.025 


3.20e-03 


7.28 


0.73 


539 


0.5 


0. 


10 


0.025 


2.40e-03 


5.46 


0.55 


534 


0.5 


0. 


10 


0.050 


4.00e-03 


9.11 


0.91 


535 


0.5 


0. 


10 


0.050 


3.20e-03 


7.28 


0.73 


536 


0.5 


0. 


10 


0.050 


2.40e-03 


5.46 


0.55 


547 


0.5 


0. 


10 


0.050 


1.80e-03 


4.10 


0.41 


546 


0.5 


0. 


10 


0.050 


1.20e-03 


2.73 


0.27 


540 


0.5 


0. 


10 


0.100 


4.00e-03 


9.11 


0.91 


541 


0.5 


0. 


10 


0.100 


3.20e-03 


7.28 


0.73 


542 


0.5 


0. 


10 


0.100 


2.40e-03 


5.46 


0.55 


549 


0.5 


0. 


10 


0.100 


1.80e-03 


4.10 


0.41 


548 


0.5 


0. 


10 


0.100 


1.20e-03 


2.73 


0.27 


614 


0.5 


0. 


10 


0.200 


4.00e-03 


9.11 


0.91 


615 


0.5 


0. 


10 


0.200 


3.20e-03 


7.28 


0.73 


616 


0.5 


0. 


10 


0.200 


2.40e-03 


5.46 


0.55 


617 


0.5 


0. 


10 


0.200 


1.80e-03 


4.10 


0.41 


696 


0.5 


0. 


10 


0.200 


1.20e-03 


2.73 


0.27 


689 


0.5 


0. 


10 


0.300 


4.00e-03 


9.11 


0.91 


688 


0.5 


0. 


10 


0.300 


3.20e-03 


7.28 


0.73 


639 


0.5 


0. 


10 


0.300 


2.40e-03 


5.46 


0.55 


698 


0.5 


0. 


10 


0.300 


1.80e-03 


4.10 


0.41 


697 


0.5 


0. 


10 


0.300 


1.20e-03 


2.73 


0.27 
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